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Abstract. We present an analysis of dust grain emission in the diffuse interstellar medium of the Small Magellanic Cloud 
(SMC). This study is motivated by the availability of 170 /jm ISOPHOT data covering a large part of the SMC, with a resolution 
enabling to disentangle the diffuse medium from the star forming regions. After data reduction and subtraction of Galactic 
foreground emission, we used the ISOPHOT data together with HiRes IRAS data and ATCA/Parkes combined Hi column 
density maps to determine dust properties for the diffuse medium. We found a far infrared emissivity per hydrogen atom 30 
times lower than the Solar Neighborhood value. The modeling of the spectral energy distribution of the dust, taking into account 
the enhanced interstellar radiation field, gives a similar conclusion for the smallest grains (PAHs and very small grains) emitting 
at shorter wavelength. Assuming Galactic dust composition in the SMC, this result implies a difference in the gas-to-dust ratio 
(GDR) 3 times larger than the difference in metallicity. This low depletion of heavy elements in dust could be specific of the 
diffuse ISM and not apply for the whole SMC dust if it results from efficient destruction of dust by supernovae explosions. 
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1. Introduction 

The Small Magellanic Cloud (SMC) is a dwarf irregular ga laxy 
at a distance of approximately 60 kpc JWesterlundlli99ll) . Its 
relative proximity enables to have a global view of the galaxy, 
as well as a sufficiently good resolution to disentangle the 
diffuse medium from the star forming regions and molecular 
clouds. The SMC is an ideal target to study the low metallicity 
diffuse medium present not only in dwarf galaxies, but also in 
the external parts of large spirals. 

The SMC contains severa l Hn regions and molecu- 
lar clouds ( Rubi o et al ] Il993al) . It has a high interstellar 
radiation field (ISRF), between 4 and 10 times that of 
the Solar neighborhood (lLeaueuxiri979t IVangioni-Flam et ID 
ll980t lA zzopardi et al. 198 8I). It is extr emely rich in neutral 
atomic gas dStaveley-Smith et alJI 19981) and has a low heavy 
element abundance jDufouii l 1984: Sa uvage & Vigrouxl ll991). 
All of these characteristics are indicative of a " young" galaxy, 
actively forming stars, which might in some sense be consid- 
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ered as a local model of distant galaxies at the beginning of 
their evolution. 

Dust grains are composed mainly of heavy atoms (e.g. car- 
bon, silicium, oxygen, . . . ). In a low metallicity environment, 
such elements are lacking. How does this deficiency act on the 
dust grain abundances? The gas-to-dust ratio (GDR) is often 
assumed to be proportional to the metallicity (i.e. there is a 
fixed fraction of heavy elements in dust). Based on a mo del of 
the ev olution of elemental abundances in gas and dust, iDwekl 
dl998l) concludes that this proportionality applies across the 
disk of the Galaxy. However, observations in dwarf irregu- 
lar and blue compact dwarf galaxies show a large scatter in 
the GDR versus metallicity relation. This indicates that a sim- 
ple scalin g might be a too simplistic description of the phe - 
nomenon lilsa et all 1 19901 Lisenfeld & Ferraralll997t [l998). 
The GDR is expected to vary among ISM components because 
destruction of dust by shock waves is more efficient in the dif- 
fuse gas while accretion of heavy elements on dust dominates 
in molecular clouds. The integrated GDR is thus expected to 
depend on the relative fraction of matter in diffuse/dense com- 
ponents and on the supernova r ate and thereby on the star for- 
mation history (Hirashi tal 19991) . 

The goal of our study is to characterize the dust properties 
in the diffuse medium of the SMC, and in particular the GDR. 
Comparing the GDR and the dust spectral energy distribution 
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for the SMC and the Galaxy, we quantify the effect of the dif- 
ference in metallicity on dust abundance and size distribution 
(e.g. the relative amounts of PAHs, very small grains and big 
grains). 

Since the gas and the dust masses are closely mixed, where 
the radiation is homogeneous the dust emission is expected 
to correlate with the distribution of the gas. In the diffuse 
medium, gas is mainly composed of neutral hydrogen and the 
GDR can be derived from the comparison between the dust 
infrared emission with the Hi column density, provided that 
the dust temperature is known. This is not the case in star- 
forming regions, where the infrared emission from dust in H n 
regions and molecular clouds closely associated with luminous 
stars can become dominant. Several authors have computed 
GDRs, comparing the IRAS infrared emission with the H i col- 
umn density, taking the 100£fm/60£fm inte nsity ratio to derive 
the dust temperature (Schwering & Israel 1989; Sauvage et al. 
ll99Cj IStanimirovic et alJkOOd but this ignores the contribu- 
tion from small grains at 60/mi. An additional difficulty comes 
from the variations in the intensity of the radiation field inside 
a beam, which casts doubt on any GDR estimate from stud- 
ies where the angular resolution is too coa rse to separate star 
forming regions fro m the diffuse emission (Li & Draind l2002: 
lAguirreetalfeoni 

The availability of 170 /mi ISOPHOT observations en- 
ables us to have a new vie w on dust pr o pertie s of the SMC. 
Independently of our work, IWilke et all d2003h have also re- 
duced and analyzed these dat a but with a distinc t focus on point 
sources. In a second paper ( Wil ke et al.ll20 04). these authors 
adress the same question as us but they investigate the main 
body (bar) of the SMC without discarding the star forming re- 
gions. 

After data reduction and validation (Sect. [2] and [3}, the 
170 //m map is compared to additional observations (Sect.|4j, in 
particular the IRAS High Resolution map (HiRes) at 100 /mi. 
Having obtained the 170/100 /mi emission ratio, in Sect. 15.11 
we derive a reference temperature for the big grains in the dif- 
fuse medium of the SMC. A n Hi column density m ap from 
ATCA/Parkes combined data ( Stanimi rovic et al.ll999t) is used 
to study the dust to gas correlation, allowing us to estimate 
the emissivity of the dust per hydrogen atom for the diffuse 
medium (Sect l5.2t . Finally, the spectral energy distribution of 
the d ust grains in the di ffuse medium of the SMC is fitted with 
the iDesert et al.l (|l990) model to quantify the abundances of 
the different grain components (Sect. l5~4l . The results are dis- 
cussed in Sect.|6l Sect.[7]presents our conclusions. 



2. The Data 

This study relies on archive data. Some were ready-made as 
the IRAS maps at 12, 25, 60, and 100 //m, or the ATCA/Parkes 
combined data of the Hi line at 21 cm. The ISOPHOT data 
was imported from the ISO Data Archive 1 . But the quality of 
the ISOPHOT map at 170 /mi in the ISO archive was not suf- 
ficient for this study, so we retrieved, reduced and calibrated 



the data. We first present the ISOPHOT data and the reduction 
processes. Then we describe the ancillary data that were used. 

2.1. The ISOPHOT data 

The Small Magellanic Cloud was almost fully covered by 
10 ob servations with the ISOPHOT inst rument jLemke et al] 
1996 D on board the ISO satellite jKessler et al] Il996: 
Laur en's et alJl2003l) . The observations were made using the 
P22 AOT (Astronomical Observation Template), i.e. a raster 
map on a two-dimensional regular grid, with the C200 detector 
(2x2 pixels) at an effective wavelength of 170 /mi. The pixel 
field of view is 1.5 ' and the raster pointing spacings are 180" 
in both in-scan and cross-scan directions. The ISOPHOT point 
spread function can be approximated by a 2D gaussian with 
90 " full width half maximum (FWHM) in both directions and 
is supposed to be constant across the map. 

We reduced the data with PIA 2 (ISO PHOT Interactive 
Analysis) V10.0 JGabriel et alJl997l fl998). except for the flat- 
field correction. In the following, we describe the specific treat- 
ments applied. Some observations showed specific problems on 
raw data (e.g. a pixel went into erratic state or a severe glitch 
affected a pixel until the end of the observation). We choose to 
remove the data for these pixels in the corresponding observa- 
tions. If this removal is not done, the reduction processes will 
interpret these erratic fluxes as flat-fielding. This distorts the 
photometry of the other pixels for the whole observation. 

To deal with the discarded pixels, we choos e to perform the 
flat-fie ld correction without the PIA tool as in Lagac h~& Dole! 
(2001). The flat-field correction is on average 1.01 + 0.05, 
0.93 + 0.03, 1.17 + 0.09 and 0.92 + 0.03 for pixel 1, 2, 
3 and 4 respectively. Reproducible flat-field correction val- 
ues for the IS OPHOT C_200 array have been published by 
Lagache &Dolel J200ll) for the carefull reduction of 3 deep 
fields (FIRBACK). We do not observe the same average flat 
field corrections factor and we find a higher dispersion from 
one observation to another. This difference may be because 
the responsivity of ISOPHOT C_200 pixels changes with the 
observed intensity, so that the flat-field correction computed 
with low fluxes is no longer valid for the brightest pixels 

J ; II 1 ° or 

(Pagani et al. 2003). This effect might affect the 170 /mi flux 
determination for bright sources. The repro ducible behavior of 
pixels observed bv lLagache & Dold J200lh is then due to the 
relative flatness of deep fields. 

The different observations have been done at different 
times, therefore the "absolute calibration" performed using the 
FCS measurement can be slightly different from one raster to 
another. This difference is observed as discontinuities between 
maps in the mosaic. To correct for these small absolute cali- 
bration differences (for a border between two maps, the mean 
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2 The ISOPHOT data presented in this paper were reduced us- 
ing PIA, which is a joint development by the ESA Astrophysics 
Division and the ISOPHOT Consortium with the collaboration of 
the Infrared Processing and Analysis Center (IPAC). Contributing 
ISOPHOT Consortium institutes are DIAS, RAL, AIP, MPIK, and 
MPIA. 
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difference in brightness observed i s about 13%), we applied a 
" general flat-field" correction as in Lagache & Dold 12001). 

After concatenating the timelines, the raster is projected on 
the sky using each individual pixel coordinates, with a pixel 
size for the map of 10 " to match the IRAS HiRes gridding. 
Because the raster step size of each observation was a full de- 
tector size, and the spacecraft Y-axis was not parallel to right 
ascensions or declinations, we observe gaps between the indi- 
vidual array pointings. The final map is presented at Fig. Rafter 
a bilinear interpolation to reconstruct the signal in the gaps. 

2.2. Ancillary data 

2.2.1. IRAS HiRes data 

As a complement to the ISOPHOT data at 170 yum, we used 
IRAS high resolution (HiRes) data at 12, 25, 60 and 100 yum. 
Maps of 2° x 2° centered at 00 A 57"'29.12* -72°32'29.2" in 
J2000 were obtained. The pixel size is 10 ". The HiRes point 
spread function (PSF) varies spatially in a map. A gaussian ap- 
proximation of the PSF is given with the HiRes maps and is de- 
fined by two spatially dependent FWHMs and a rotat ion angle. 
For m ore in formation on HiRes, pl ease refer to lAumann et alJ 
(1990) and Bontekoe et al. (11994 . Since we are searching 
for extended structure propert ies, the H iRes IRAS maps were 
brought to the DC-calibration JWheelock et alJl994l) . This DC 
calibration is only valid for spatial scales greater than 40 '. As 
the SMC exhibits structure on various spatial scales, the DC 
calibration introduces some limitations that will be discussed 
in Sect.l3~2l 

2.2.2. ATCA/Parkes combined SMC H i data 

To study the correlation between gas and dust in the SMC, 
we needed a Hi column density map to compare with the 
ISOPHOT map, in order to compute the abundances of dust 
gra ins relative to the hydroge n . 

IStanimirovic et all dl999l) have combined Parkes telescope 
observations of neutral hydrogen in the SMC with the Australia 
Telescope Compact Array (ATCA) aperture synthesis mosaic 
to obtain a set of images sensitive to all angular spatial scales 
between 98 " and 4° from 90 to 215 km.s -1 . The spatial resolu- 
tion obtained is 98 " and the velocity resolution is 1.65 km.s -1 . 
We used a column density map of the SMC, obtained by in- 
tegrating over the whole velocity range and assuming that the 
neutral hydrogen is optically thin 3 . 

2.2.3. Parkes H i data for the Galactic foreground 

A Galactic filament exists in front of the SMC. Even though 
it is very smooth, large (see contours of Fig. [0 and uniform 
(column densities between 2.5 • 10 and 4 • 10 atom.cm ), 
it can bias studies. The far-infrared fluxes observed in either 
the ISOPHOT or the IRAS HiRes maps represent the emission 
from the SMC dust, but also from this Galactic foreground. 

3 The column density map and the data cube are available at 
ftp://ftp2.naic.edu/pub/ast/sstanimi/. 



However, this Galactic component of the emission is not ob- 
served in the ATCA/Parkes data, due to a velocity cut and spa- 
tial filtering. To study the SMC dust properties, it was thus nec- 
essary to remove this foreground component (see Sect. I3.3i . 
In this step, we used a Galactic H i column density map con- 
structed from data of the Parkes narro w-band H i survey of the 
tidal arms of the Magellanic system (Bri ms et alJl20 00). The 
column density map of the Hi Galactic foreground was ob- 
tained by integrating the emission from -60 to +50 km/s, out- 
side the SMC velocities. The resolution is 15'. Due to side- 
lobes effects, the map obtained can slightly overestimate the 
column density by 10 to 20 %. 

3. Photometric consistency 

DIRBE provided infrared absolute sky brightness maps in 10 
bands in the wavelength range 1.25 to 240 fim. This abso- 
lute photometry is used to check ISOPHOT mosaic and IRAS 
fluxes. 

3.1. ISOPHOT data 

We checked if the ISOPHOT 170 /mi flux is consistent with 
the spectral energy distributions obtained from the 100, 140 
and 240 /im DIRBE fluxes in the SMC. To do this we con- 
volved the ISOPHOT map with the DIRBE beam. This gives 
18 ISOPHOT fluxes in the SMC as if they were observed by 
DIRBE and we compare ISOPHOT and DIRBE fluxes by plot- 
ting the spectral energy distribution and fitting it with a modi- 
fied black body with a spectral index of 2 on the DIRBE val- 
ues. Fig. shows that the DIRBE-convolved ISOPHOT flux 
integrated over the 18 positions observed by DIRBE, is con- 
sistent within the errors bars with the best modified black- 
body fit (obtained for a temperature of 20. 6K, an emissivity 
of e#(250/um) = 3 • 10~ 28 cm 2 for 10 20 hydrogen atoms. cm -2 ). 
Please note that these values correspond to an average over the 
whole galaxy and were obtained only for calibration purposes. 

We have also plotted the individual fit values at 170 /vm 
versus the corresponding 18 ISOPHOT values. This plot shows 
a significant scatter but is consistent with a linear correlation 
with unity slope and zero offset. 

A comparison with the IRAS high resolution data was also 
performed to check possible photometric errors introduced by 
the interpolation process. The correlation between IRAS and 
ISOPHOT data does not change when we restrict the compar- 
ison to the pixels measured by ISOPHOT rather than the full 
interpolated map. The interpolation required to fill the holes 
in the measurements does not seem to bias the image compar- 
isons. 

3.2. IRAS data 

The IRAS HiRes maps have the same photometric calibration 
as the IRAS ISSA maps. The IRAS detector had a dwell-time 
dependent responsivity change: the gain changes as a function 
of source size. This effect is known as the " AC/DC effect". 
A first set of corrections were taken into account in the IRAS 
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Fig. 1. The ISOPHOT map at 170 fim (in MJy/sr) after reduction processes and a bilinear interpolation to fill the holes left by 
the mapping pattern. The flat-field has been efficiently corrected for the diffuse part of the map, but one can still see low level 
defects in the brightest parts. Contours of the Hi Galactic foreground are overplotted at column densities of 2.6, 2.8, 3., 3.2, 3.4 
and 3.6- 10 20 H atom.cm-2 (from light grey to black). 



data products. These corrections were later found to be insuf- 
ficient. On the basis of a comparison of the IR AS and D IRBE 
sky brightness for extended emission, IWheelocketai] i ll 9941) 
recommended applying additional correction factors of 0.88, 
1.01, 0.82 and 0.74 at 12, 25, 60 and 100 fim respectively. 

However, these correction factors are only valid for spatial 
scales greater than 40'. This is not true in the S MC where struc- 
tures o n various spatial scales are present. Stanimirovic et al. 
(2000) have also compared DIRBE and IRAS data by inte- 
grating fluxes over 6.25° 2 in the SMC. They found that inte- 
grated flux densities measured by DIRBE for 60 and 100 fim 
are higher than IRAS fluxes by about 10 to 20%, which is a n 
opposite correction to that observed bv lWheelock et al] ([1994). 

We have compared the IRAS (convolved to the DIRBE res- 
olution) and the DIRBE 100 fim maps on a 10° x 10° scale 
arou nd the SMC to und erstand this discrepancy. We see that the 
IWheelocketai] ([l994) correction factor of 0.74 is only truly 
valid outside the galaxy and increases in the SMC up to a max- 
imum value in the north of the o ptical bar (N66 region) . For the 
diffuse medium of the SMC, the Whe elock et alJlll994l) correc- 
tion factor is sufficient. After correction, we estimate that the 



100 fim IRAS brightness for the diffuse medium is underesti- 
mated by less than 15%. 

3.3. Galactic foreground removal 

Unlike the ATCA/Parkes Hi map, the far-infrared emission 
maps (ISOPHOT and IRAS) trace both the dust emission in the 
SMC and the dust in our Galaxy. To disentangle these two com- 
ponents in the far-infrared emission, we used a column density 
map of the Hi Galactic foreground presented in Sect. 12.2.31 

Dust grains absorb energy from the UV/optical radiation 
field that they re-emit in the infrared wavelength as thermal 
emission. Large grains are at a fixed temperature and their emit- 
ted intensity can be expressed as: 



N H e H {v)B v {T dust ) 



(1) 



where Nh is the column density of hydrogen and e//(v) is the 
emissivity of grains per hydrogen atom. Assuming that the 
emissivity of dust grains in the far-infrared follows a power 
law: 



e H (v) = e H (v Q ){-f 
vo 



(2) 
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Fig. 2. Comparison of the DIRBE-convolved ISOPHOT fluxes 
with the DIRBE ones (in MJy/sr), after an integration on the 
whole SMC. The modified black body traced was fit on the 
DIRBE values only and correspond to a temperature of 20.6 K 



and an emissivity of e//(250/mi) = 3-10 28 cm 2 for 10 2U hydro 
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gen atoms. cm . We see that the ISOPHOT value at 170 /mi is 
consistent with the DIRBE photometry. 



and taking (3 = 2 (as found for Galactic emission 
jLagache et aljfl998l 1 199 9) for wavelengths between 100 /mi 
and 500 /mi). 

Using a temperature of 17.5 K and an emissivit y per hy- 
drogen atom of 10~ 25 cm 2 (Boulanger et alJll996al) 4 for the 
Galactic dust, we computed the Galactic thermal emission 
and removed this foreground emission from the IRAS and 
ISOPHOT maps to obtain the emission of the SMC dust only. 
At 12, 25 and 60 //m, we used the Inlhao, hs/hoo an d I&alhw 
ratios computed from IRAS maps outside the SMC, together 
with the computed foreground emission map at 100 /zm to re- 
move the Galactic foreground emission at these wavelengths. 
However, this foreground emission remains in the data noise at 
12, 25 and 60 /mi. 



4. Spatial image comparison 

4.1. Projection and resolution 

The data set presented in Sect. |2] is composed of maps with 
different fields of view and different resolutions. We projected 
all of them onto a common grid: that of IRAS HiRes maps 
(2° x 2° field of view centered at 00' 1 57 m 29.12 s - 72°32'29.2" 
(J2000) with 10 " pixel size). These maps were inter-compared 
by pairs. 

To compare the ISOPHOT map with one of the HiRes 
maps, the projected and re-sampled ISOPHOT map was con- 
volved with a spatial varying kernel in order to be at the HiRes 
resolution when the point sprea d function of I SOPHOT was 
smaller than the one of HiRes jKerton & Martini 1200 lb . The 
HiRes map was also convolved with a spatially varying kernel 
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to match the ISOPHOT resolution when the HiRes beam was 
smaller than the ISOPHOT one. 

To compare the HI column density map with the ISOPHOT 
map, no convolution was necessary as long as their resolutions 
were comparable. 

Fig- E shows two of the transformed maps: the ISOPHOT 
map degraded to the IRAS 100 /mi resolution and the projected 
ATCA/Parkes map. 

4.2. Colors from scatterplots 

To characterize the diffuse medium (i.e. the extended region 
where the infrared emission-Hi column density correlation 
is homogenous) spectral energy distribution, we developed a 
comparison method that is not biased by point sources and is 
also not affected by offsets in the data calibration. 

We work with two maps (obtained after the treatment pre- 
sented in Sect. @}. The general correlation between the two 
maps is first quantified by fitting a straight line to the scatterplot 
of the pixel-to-pixel comparison. We deduce a first ratio value 
and an error cr given by the dispersion in the correlation, which 
enables us to filter out localized sources (i.e. star-forming re- 
gions) by removing pixels that are 3<x above the correlation 
(there is no points below — 3cr). 

From the " filtered" maps, we compute a ratio map by lin- 
early fitting the pixel correlation in the scatterplot for a sliding 
box of 50x50 pixels. This corresponds to a 5.5x5.5 beams box 
which is a minimum size to have independent data points. The 
size of the box could not be greater, or we would not trace the 
local variations. In some cases there are too few points inside a 
box (too many points above 3cr) to do a meaningful fit. We as- 
sign no value to these positions. We also skip positions where 
the correlation coefficient is too low. Adjusting a gaussian on 
the histogram of the ratio values gives us a reference ratio and 
the associated dispersion. This dispersion contains photomet- 
ric and fitting errors but also a true scatter in the ratio. Both 
phenomena cannot be easily disentangled. 

The spatial distribution and the histogram with the gaussian 
fit obtained for the /ioo/im/ hiO/im ratio are presented in Fig.|4]as 
an example. We see that the comparison method avoids large 
areas in the main body (bar) of the SMC. These regions corre- 
spond to highly structured star-forming sites. In these regions, 
either the filtering technique has left too few independent posi- 
tions in one box to compute a significant correlation or no clear 
correlation is observed. 



5. Results 

5.1. Infrared colors 

The method presented in Sect. I4.2I was applied to the IRAS 
HiRes and ISOPHOT maps to compute the /ioo/im/^i70^m, 
liatimlhiQiim, hstanlhiQism, hipim/ hio^m ratios. Adjusting a 
gaussian on the histograms gives us reference ratios and their 
associated dispersions: 



We checked that the observed IRAS and DIRBE far-infrared ratios 



around the SMC are consistent with these values. 



'100/itn 
hlOfaa 



= 1.0 + 0.2 



(3) 
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Fig. 3. ISOPHOT map (left panel, in MJy/sr) and ATCA/Parkes H i column density map (right panel, in 10 20 Hydrogen atom/cm 2 ) 
after projection on the IRAS HiRes grid, and at the same resolutions. For the ISOPHOT map, logarithmic contours are overplotted 
at fluxes 11.2, 14.1, 17.8, 22.4, 28.2, 35.5, 44.7, 56.2, 70.8 and 89.1 MJy/sr. 



hlOfim 



0.5 + 0.1 



0.02 + 0.01 



= 0.011+0.006 



(4) 
(5) 
(6) 



Using equation^ and [2] we can derive from the rat i 

'170 /im 

map the big grain equilibrium temperature distribution in the 
diffuse medium in the Small Magellanic Cloud. This leads to a 
reference temperature of: 7d us t = 22 ± 2 K. Note that this tem- 
perature is different from the one that would be i nferred f rom 
the r atios of the mean 7ioo /zm and /170 i_, m fluxes. IWilke etafl 
(2004) found a mean temperature for the main body of the SMC 
of 20.9 ± 1 .9K, which is lower than the temperature we find for 
the diffuse medium. Their mean value falls within our range of 
temperatures, which is expected because the regions we study 
overlap. The difference in mean values could be due to the pres- 
ence of cold atomic and molecular clouds in the main body of 
the SMC. Our studies seem to be consistent with each other and 
are more complementary than comparable. 

The large dust grain temperature is related to the radiation 
field intensity and will be discussed in Sect. [6] 

5.2. Gas to dust correlation 

From the ISOPHOT and the ATCA/Parkes maps, using the 
method described Sect. 14. 21 we found the distribution presented 
on Fig. [5] This distribution exhibits a main peak and a distribu- 
tion tail. The peak is characterized by: 



200 - 




o.s 



'170//m 



in 



= 0.17 • 10~ 20 ± 0.06 ■ 10~ 20 MJy.srKcm 2 



(7) 



Fig. 5. Histogram of the 170 micron flux to Hi column den- 
sity ratio (in MJy.cm 2 /sr). The gaussian fit and the residuals are 
presented as continuous and dotted lines respectively 



Since the regions used to define the ratios are not exactly 
the same, the determination of the characteristics of the diffuse 
medium could be biased. We checked that this is not the case. If 
we use the region where Inofim/Nm is computed to determine 
the colors, we find that the derived quantities are not affected 
by more than 15%. We see that the characteristics of the dif- 
fuse medium are not biased by the definition of the limits of 
the studied region. Comparing the region defined by equation 
[7] with the Hd distribution, we see in Fig. |6] that the regions 
corresponding to the tail in the Iho^/Nhi (located mainly in 
the main body) are regions where the diffuse H a emission is 
important. Comparing the map in Fig.|4]and Fig. [6] we also see 
that the regions where no ratios can be computed correspond to 
important star-forming sites. 
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Fig. 4. Histogram and spatial distribution of the 100-170 fim colors. For the histogram, the gaussian fit and the residuals are 
presented as continuous and dotted lines respectively. For the map of the 100-170 fim colors, contours correspond to the 0.8 
(black) and 1.2 (white) levels (one sigma limits of the gaussian fit). It should be noticed that the noise in the color map is not 
homogeneous. 




Fig. 6. Uncalibrated Ha map. As an illustration, the contour 
line corresponds to /no^m/A^/ = 0.23 MJy.cm 2 /sr (see equa- 
tion0 as in figure^] 



Assuming that all the dust is in the Hi, we use the 
hiQiimlNm ratio map and the temperature map computed above 
to obtain an histogram of the emissivity per hydrogen atom of 
the dust grains in the diffuse medium. We find a reference emis- 
sivity at 170 //m of: 



e H (A) = 7.3 ■ 10~ 27 + 3.2 • 10^'cm 



-27 „ 



(8) 



The error bar takes into account uncertainties in the tem- 
perature and /no fmx/Nffi- If the emissivity follows a power 
law with a spectral index of 2, the emissivity per Hi 
atom in the SMC at 250 ^m is e H (A) = 3.4 • 10~ 27 + 



1.5 • 10~ 27 (/l/250 fimy 2 cm 2 . This corresponds to k v — 
0.21(/l/250/iirfr 2 cm 2 g~'. Assuming similar dust optical prop- 
erties in the SMC and the Galaxy, the emissivity is related to 
the GDR and the ratio between hydrogen column density and 
E(B - V) color excess by: 



e H {SMC) 
eniGalaxy) 



M^l(SMC) 



N(H,)+N(H 2 ) 
E(B-V) 



{Galaxy) 



— (Galaxy) 



N(Hj)+N(H 2 ) 
E(B-V) 



(SMC) 



(9) 



The solar neighb orhood value for the emissivity is 
10~ 25 cm 2 at 250 jum feoulanger et al.lll996ah . The inferred 
SM C GDR is th u s 30 times greater than in the Galaxy. 
IStanimirovic et all d2000l) have presented a similar study us- 
ing the same Hi data together with HIRAS data. They came 
to the same conclusion that the GDR is a factor 30 higher in 
the SMC than in our Galaxy. However, this agreement might 
be a coincidence since 1) we did not use the same calibra- 
tion of the IRAS data, 2) their temperature estimation based 
on the 60/100/im color is higher than ours, 3) the GDR they 
present is averaged over the whole galaxy (including star form- 
ing regions) and 4) their GDR is compared to a Galactic GDR 
of 375 instead of the standard value of 100 5 . Using the stan- 
dard Galactic GDR valu e, the average SMC GDR found by 
IStanimirovic et al]j2000h becomes 100 times higher than in the 
Galaxy. 

Our estimate of the SMC's GDR 30 times higher than 

in the Galaxy, is significantly lar ger than the d ifference in 

metallicity, which is less than 10 (Dufoug ll984l) . It is twice 

the observed difference in the N ^}„ N S^ 2) values presented by 

1 J ^3. t(b-v) r J 

Tumlinson et al. (2002). 



This value can be ch ecked with the interstel lar a bundance stan- 
dards ta ken from table 2 of lSofia & MeveJ feOOlbh and lSofia & Meved 
1 200 3). 
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Fig. 7. Spatial distribution of the hio umlNni ratio (in 
MJy.cm 2 /sr). Regions where no ratio is computed (see Sect. 
14.21 are left blank. The contour line corresponds to a 
hio umlNm = 0.23 MJy.cm 2 /sr level. We see that the IR ex- 
cesses are concentrated near the optical bar. 
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Fig. 8. Best fit for the dust infrared emission spectrum, ob- 
tained with the Desert et al. model. The brightness corresponds 
to an interstellar gas column density of 10 20 cm~ 2 . The IRAS 
and ISOPHOT data are represented by squares. The fit obtained 
with the model (continuous line) corresponds to the emission 
of three components: PAHs(dashed line), very small grains 
(mixed line) and big grains (dot-dot-dot-dashed line). 



5.3. Infrared excesses 

A noticeable result of the ISOPHOT and H i comparison is the 
presence of a tail in the 170 fxm/H i histogram (see Fig. [3}. This 
tail corresponds to pixels in the SMC optical bar with excess 
170 fim emission with respect to Nhi (see Fig. [7}. We checked 
that the mean temperature for this region of excess is the same 
as that for the diffuse medium. These excesses may be due ei- 
ther to a variation in the dust abundances between the diffuse 
medium and the bar, or due to the presence of cold H i or H2 
gas. We believe that the excesses are real but it is difficult to 
infer an amount of cold H 1 and molecular gas from it. Firstly, 
it is difficult to disentangle these two possible contributions. 
Secondly, where some of these excesses are due to dust in cold 
H 1 or H2 we do not expect them to be correlated with H 1 emis- 
sion and another method is needed to quantify them. It is thus 
beyond the scope of this paper to translate these excesses into 
an estimate of cold gas that could be compared to other stud- 
ies where this is attempted ( Stanimiro vic et alJ l2000). Finally, 
we feel that the photometric uncertainties in the ISOPHOT 
and IRAS data are an additional difficulty, especially near star- 
forming regions. Further studies with the Spitzer satellite could 
address this question with better photometry and angular reso- 
lution. 

5.4. SMC dust properties in diffuse H 1: modeling the 
spectral energy distribution 

The different intensity ratios obtained in Sect. 15.11 enable us 
to trace the spectral energy distribution of dust in the diffuse 
medium of the SMC and to compare it with that of the so- 
lar neighborhood. The difference in temperature and in abun- 
dances is seen Fig.|9]as a shift along the wavelength and bright- 
ness axis respectively. The error bars include true variations in 



the dust properties across the SMC diffuse interstellar medium 
which are correlated over wavelengths. 

We used the dust model of Desert et al. (1990) to quan- 
tify the grain dust abundances. In order to model the spectral 
energy distribution (SED), we need to estimate the interstellar 
radiation field (ISRF) in the diffuse medium of the SMC. In 
Sect. 15. II we found the big dust grain equilibrium temperature 
to be 22 K. Using the dust thermal equilibrium temperature in 
the solar neighborhood of 17.5 K (Boulang er et alll996al) as a 
reference, the dust temperature in the SMC leads to a radiation 
field xsmc = T 4 d ^,(SMC)/T 4 d ^ t (Galaxy) = 4. This value is 
consistent with the estimate of t he IS RF made from massive 
star counts by iLeaueuxl dl979l) . IVangioni-Flam et a 
also found the same ratio from UV measurements. We used the 
Mat his et aD (1983) 10 kpc Galactic ISRF and scaled it by a 
factory = 4 for the modeling. We have adjusted only the abun- 
dances of PAHs, very small grains (VSG) and big grains. For 
the other parameters defining the size distribution, we have kept 
the Galactic values. The best fit is presented on Fig. [8] it corre- 
sponds to under-anbundant PAHs, very small grains (VSG) and 
big grains (BG) compared to the Galaxy (36, 23 and 30 times 
less respectively). Mass abundances relative to hydrogen, com- 
pared with Galactic values (Desert et al. 1990) are summarized 
in table [T] 

The modeling confirms the GDR estimated in Sect l5.2l to be 
around 30 times greater than in the Galaxy. Within the uncer- 
tainties, PAHs and very small grain abundances are also con- 
sistent with this factor. 

The model is clearly below the 60 /itn point. This could 
be due to our crude estimate of the ISRF. We have used a sin- 
gle ISRF while our data analysis shows temperature variations 
across the SMC diffuse medium. To take the variations in the 
radiation field into account, we transformed the temperature 
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Fig. 9. Infrared spectral energy distribution of dust in the diffuse medium for the SMC and the Solar Neighborhood. The bright- 
ness corresponds to an interst ellar gas column den sity of 10 20 cirT 2 . The Galactic spectrum (in grey) is composed in the far-IR by 
COBE/FIRAS data from Boulanger et al. ( 1996a) (with the corresponding spectrum from grains in thermal equilibrium indicated 
by the dotted line), and at shorter wavelengths by the COBE/DIRBE (Bernard et al. 1994) and Arome balloon data (triangle at 
3.3 fj.m) dGiardet al.ll 1994 . An ISOCAM spectrum dBoulanger et alJll996bl) scaled to the 12 fim brightness is also shown. The 
SMC spectrum (in black) is the best fit obtained with the Desert et al. model, same as in figure[8] 



Table 1. Mass abundances relative to hydrogen in the SMC 
and the Galaxy. The grain size distribution is assumed to be the 
same in both cases. 



Component 


SMC 


Galaxy 




m„ 




m„ 




PAH 


1.2 


io- 5 


4.3 


1(T 4 


VSG 


2.0 


1(T 5 


4.7 


1(T 4 


BG 


2.2 


icr 4 


6.4 


io- 3 



distribution obtained in Sect. 15. ll into a. xsmc distribution. For 
each value of the distribution, we computed a sp ectral energy 
distri bution (SED) of the dust emission with the D esert etall 
(1990) model, using the Mathis 10 kpc ISRF multiplied by the 
corresponding^ as an input Of is a constant for the whole spec- 
trum). We then computed the resulting total SED by adding the 
individual SEDs with weights according to the x distribution 
shape. The resulting fit is not sufficiently broadened to encom- 



pass the 60 /jm point. In addition, the mass abundances of dust 
grains have to be reduced by 8, 5 and 7 for the PAH, VSG and 
BG components respectively. Due to the high non-linearity of 
the Planck curve, even a small broadening of the x distribu- 
tion tends to reduce the abundances of dust grains significantly. 
These values are thus indicatives only but they show that the 
mass abundances represented in tableware more likely to be 
upper limits and that the emissivity could be even lower. 

It is possible to fit the 60 /urn value with a single ISRF scaled 
by X - 10 (see Fig. 1101. In this case, the dust grain mass abun- 
dances must then be reduced by 3, 2 and 2 for the PAHs, very 
small grains and big grains respectively. Such a high ISRF cor- 
responds to a dust temperature of 25.6 K, which is higher than 
the estimate obtained in Sect. l5.1l but still compatible with the 
100 /urn and the 170 /urn fluxes. This higher temperature would 
lead to an emissivity of 2.4 • 10~ 27 cm 2 at 250 /urn, which is even 
lower than the estimation obtained in Sect. 15.21 and results in 
a GDR 40 times higher than in the Galaxy. A higher radiation 
field could thus account for the 60yum excess without contra- 
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Fig. 10. Dust infrared emission spectrum fit obtained with the 
Desert et al. model to adjust the 60/im excess. The IRAS and 
ISOPHOT data are represented by squares. This spectrum is 
obtained with a single ISRF scaled by \ — 10 an d reduced 
dust grain mass abundances for PAHs (dashed line), very small 
grains (mixed line) and big grains (dot-dot-dot-dashed line). 

dieting the main conclusion of this study, namely the fact that 
the GDR in the diffuse SMC ISM is much lower than the solar 
neighbourhood value scaled by the difference in metallicity. 

The observed 60 //m excess could also be due to a change 
in the grain size distribution (Galli ano et aljE003l) . We as- 
sumed the Galactic grain sizes distribution to apply in the 
SMC. However, by changing the relative mass contributions of 
the different grain components, we enlarge the discontinuities 
in the size distribution between the three components of the 
model. It is necessary to further increase these discontinuities 
to account for the 60 //m excess. 

Further investigations are needed to characterize the SMC 
ISRF and especially the 60 yum excess. The intensity and spec- 
tral shape of the ISRF are likely to be different than the Galactic 
one. 

6. Discussion 

Studying the gas to dust correlation in the diffuse medium of 
the SMC (see Sect. 15. 2> . we found a GDR in mass around 30 
times greater than in the solar neighborhood, larger by a fac- 
tor of 3 than the differe nce in metallicity. T his result confirms 
the value found by Stanimirovic et al. ( 2000) in a similar study. 
This lower depleti on of heavy elemen ts in dust was indepen- 
dently reported by Wel tv et alJ (|2001), who found Mg and Si 
to be essentially undepleted in th e SMC. It was also reporte d 
for a sample of dwarf irregulars bv lLisenfeld & FerradJl998l) . 
Our result in the SMC fits with the relation that they observed 
between GDR and metallicity. 

An interpretation of this relation has been proposed by 
iHirashita et all d2002l) . With a chemical evolution model, they 
illustrated the dependence of the GDR on the star formation 
history. Stars generate dust but also contribute to their destruc- 
tion by supernova shocks. They propose a scenario in which 
an intermittent star formation history leads to a time variable 



GDR. At some epochs of the galaxy evolution, the formation 
of dust grains is less efficient than the destruction process by 
supernova shocks, so that the GDR no longer scales with the 
metallicity. This model could be appli ed to the SMC w ith its 
specific star formation history (Harris & Zaritskvl l2003l) . It is 
also necessary to account for a possible difference in the GDR 
between the diffuse and the denser medium. 

Extinction studies lead to significantly smaller value of the 
GDR, closer to that expected fro m the Galactic value for a 
linear dependence with metallicity. Bouch eTet alJ fl985) used 
photometric observations of O-B stars in the SMC to derive 
their extinction curves. They found a GDR value that is about 8 
times the Galactic value, whi ch is consistent with the difference 
in metallicity. More recentlv. lTumlinson et alJ (12002) measured 
iim^v) f° r a sam Pl e °f stars observed with FUSE and found 
a mean value 16 times lower than the solar Neighborhood 
value. This is closer but still smaller by a factor 2 than our 
estimate. The difference may be due to the fact that extinction 
studies are based on line-of-sights with high hydrogen column 
density where dense interstellar medium is likely to be present. 
Most of the stars are in the optical bar and/or in star formation 
regions associated with molecular clouds. The difference in the 
GDRs inferred in our work and extinction studies might thus 
point to a dependence on the environment: a decrease of the 
ratio from the diffuse ISM to star forming regions that might 
arise from the condensation of heavy elements on dust grains in 
dense gas and their dest ruction by supernova shocks in the dif- 
fuse ISM. Stanimirovic et al. (2000) have also noted this possi- 
bility for the SMC, but it still requires confirmation by further 
studies. 

7. Summary and conclusions 

This study was made feasible by the availability of 10 
ISOPHOT observations at 170 pun, covering most of the SMC, 
in the ISO Data Archive. After reduction processes, mapping 
and Galactic foreground removal, the ISOPHOT data were 
compared with IRAS maps and a ATCA/Parkes combined H i 
column density map to assess the properties of the dust in the 
diffuse medium of this low metallicity galaxy. 

For the dust in the diffuse medium of the SMC, we found a 
reference equilibrium temperature of 22 + 2 K and a reference 
dust emissivity of 3.4 • 10~ 27 ± 1.5 • 10~ 27 ( 2S0 ^ m )cm 2 per hy- 
drogen atom. This leads to a gas-to dust ratio 30 times greater 
than in the Galaxy. This high value reflects the SMC low metal- 
licity, but also requires a lower depletion of dust elements than 
in the Galaxy. This shows that a simple linear dependence be- 
tween the GDR and the metallicity does not apply for the dif- 
fuse medium of the SMC. 

The spectral e nergy distribution is modeled with the 
besertet al.l fl990) model. The best fit obtained shows a simi- 
lar decrease for PAHs and very small grains as for big grains. 
However, a 60 yum excess remains that can not be accounted for 
with simple assumptions. Further constraints, in particular on 
the ISRF, are needed to understand it. 

Our GDR estimate for the diffuse medium of the SMC is 
a factor 2 to 3 higher than that derived from extinction studies 
towards bright UV stars with high foreground hydrogen col- 
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umn densities. This difference supports the notion of variations 
in the GDR, from diffuse regions where grains are more fre- 
quently destroyed by supernovae, to denser and more quiescent 
regions where grains re-accrete heavy elements. 
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